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The Anderson prescription for surfaces and impurities 
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We test the Anderson prescription Q, a BCS formalism for describing superconductivity in in- 
homogeneous systems, and compare results with those obtained from the Bogoliubov-de Gennes 
formalism, using the attractive Hubbard model with surfaces and nonmagnetic impurities. The 
Anderson approach captures the essential features of the spatial variation of the gap parameter and 
electron density around a surface or an impurity over a wide range of parameters. It breaks down, 
however, in the strong-coupling regime for a weak impurity potential. 



o 

m 

"a 
o 

o _ 

u • 

Of 

3 ; 

cfl ; 



i 

-a 
c 

o 
o 



(N 
> 

o 

OS 

o 
m 
o 
o 
o 

-I— > 



i 

C 

o 
o 



X 



In microscopic treatments of inhomogeneity effects in 
superconductors, impurities are often averaged over in 
some manner In the last decade, partly because 

of the increased computational resources now available, 
and partly because of the technical advances which allow 
small particle fabrication and single atom manipulation 
H , the role of inhomogeneous effects in superconductors 
has received more wide-spread attention j| . One of the 
theoretical frameworks for addressing these questions is 
the Bogoliubov-de Gennes (BdG) formalism ||. This 
formalism allows one to answer questions regarding sur- 
faces, interfaces, and impurity effects at a level of detail 
not previously addressed. At present, however, the scope 
of problems for which one can compute results accurately 
is limited by computer resources, since matrices whose di- 
mension grows with system size require full diagonaliza- 
tion and should be solved self-consistently. On the other 
hand, the Anderson prescription first presented to 
examine the impact of impurities on superconductivity, 
is a BCS formalism for an inhomogeneous system and re- 
quires one diagonalization of the single-particle problem. 
The purpose of this paper is to examine the limits of ap- 
plicability of the Anderson approach, as compared to the 
BdG formalism. We first summarize the two approaches, 
and then follow with some concrete examples, utilizing 
surfaces and impurities as sources of inhomogeneities. 

We find that the Anderson prescription works very well 
for surfaces and in some regimes for a single impurity. 
It breaks down for strong coupling with weak impurity 
scattering. 

For our purposes we adopt a convenient model to de- 
scribe s-wave superconductivity, the attractive Hubbard 
Hamiltonian (the results of our study will presumably 
apply to d-wave or other symmetry states): 
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Here, a\ a (aj CT ) creates (annihilates) an electron with spin 
a at site i and is the number operator for an elec- 
tron with spin a at site i. The t$ is the hopping rate of 



electrons from one site to a neighbouring site (often near- 
est neighbours only are included, and we will adopt this 
model here), and \U\ is the attractive coupling strength 
between electrons on the same site. As usual, this attrac- 
tion is justified in terms of an electron-phonon coupling, 
where retardation effects are unimportant, as is the case 
with many conventional superconductors. The second 
term includes the chemical potential /i and the impurity 
potential at site i, e,. We assume that impurities act 
to raise or lower single site energy levels. It is worth 
noting that impurity effects can certainly enter in other 
ways. For example if an impurity occupies one of the 
sites, undoubtedly the hopping amplitude to and from 
that site will also be altered, as will the interaction be- 
tween two electrons occupying the same orbital on that 
site. In many studies (see, for example, Ref. 0), the 
are randomly distributed with some probability distribu- 
tion, and then the results are averaged, to reflect the fact 
that we generally have no control over the precise dis- 
tribution of impurities in the bulk. However, in systems 
where a single impurity can be added to the surface, for 
example (see, eg., Ref. 0), we would want to study this 
model with only one impurity, at a particular site (in this 
case, on the surface). 

Equation (Q) also allows us the freedom to choose peri- 
odic boundary conditions (PBC) (to recover well-known 
results) or 'open' boundary conditions (OBC). The lat- 
ter are natural in a tight-binding context; they require 
no assumptions about the order parameter, for example. 
'Open' here simply means that electrons cannot hop be- 
yond the surface. Here again more sophisticated bound- 
ary effects could be included — for example, in a real 
system the hopping integral at the surface will no doubt 
differ from that in the bulk, but we leave aside these finer 
points. 

The BdG equations are obtained by defining an effec- 
tive Hamiltonian, with effective potentials 0. By diag- 
onalizing this effective Hamiltonian through the gener- 
alized Bogoliubov-Valatin transformation || , one arrives 
at the two BdG equations B: 

E n u n (i) = ^ Mi'Unii') + V l u n (i) + AiV n (i) (2) 
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E n v n (i) = -};A u/ y n (i') - ViV n (i) + A*u n (i) (3) 

i' 

where 



A; 



Su> M - e» ■ (4) 



The self-consistent potentials, V,, and Aj, are given by 
A, = |[/| -2/„) (5) 



Vi 



|w„W| 2 /n + |«nW| 2 (l-/«) 



(6) 



We use the index n to label the eigenvalues (there are 
2N of them), the index i to label the sites (1 through 

N) , and the composite eigenvector is given by f "™ ^ , of 

total length 2N. The sums in Eqs. (|]j|) are over positive 
eigenvalues only. The /„ is the Fermi function, with ar- 
gument f3E n , where (3 = 1/ksT, with T the temperature. 
The single site electron density, rii, is given, through Eq. 
d), by Vi = -\U\nj2, 

The equations (|||]) for the effective potentials were 
determined through a variational principle so that the 
effective Hamiltonian allows fluctuations in any number 
of mean fields ||. As written, there are two possible 
mean field potentials, the Hartree potential, Vi, and the 
pair potential, A,;, from which the ground state energy 
and other properties may be obtained. 

An alternate prescription was originally proposed by 
Anderson |jj , whereby one first solves for the eigenvalues 
and eigenstates of the 'non- interacting' problem, i.e., 



E°w n (i) = y^Ag/Wnii'). 



(7) 



Using the unitary matrix, Ui n , for a basis which diagonal- 
izes the single-particle Hamiltonian, one can determine 
the transformed electron-electron interaction: 



(8) 



which now mediates the (generally off-diagonal) electron- 
electron interaction. The gap and number equations are 
obtained just as in BCS theory, except that now the la- 
bel is not the wave vector k, but rather some quantum 
number n, which simply enumerates the single particle 
eigenvalues ||. From the solution, one can obtain the 
ground state energy and, by transforming back to space 
coordinates, site-dependent quantities. 

The simplest origin of gap inhomogeneity in a super- 
conductor is the surface. The presence of surfaces beyond 
which electrons are unable to move yields a gap param- 
eter (i.e., pair potential) which can exhibit a variety of 



behaviour near the surface. Traditionally in Ginzburg- 
Landau treatments the gap function is given a priori a 
boundary condition || ; here the behaviour near a bound- 
ary (or impurity) is a derived quantity, i.e. as the solution 
to the BdG (or Anderson) equations. 

The advantage (for the Anderson approach) of exam- 
ining the impact of surfaces on the gap parameter is that 
an analytical solution exists for a simple tight-binding 



model 10 1. The eigenstates for a chain [jrT[ of length N, 
with lattice spacing a and nearest-neighbour hopping t, 
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and the eigenenergies are 
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where n = 1,2,..., A. We use these analytical results in 
the Anderson approach (making it not significantly more 
difficult than BCS theory), while in the BdG approach 
these analytical solutions are not particularly helpful. 



U|/t=1.5 N = 64 OBC 
(a) BdG 





FIG. 1. Gap parameter Ai as a function of site number i 
and electron density n, for iV = 64 with OBC and \U\ = 1.5 1. 
The BdG and Anderson results are shown in (a) and (b), 
respectively. 

In Fig. 1 we show the gap parameter A^ as a func- 
tion of site number i, for electron density n ranging from 
half-filling to zero. The chain length is 64 sites, and we 
have used OBC and the coupling strength \U\ = 1.5 1. 
Here A^ is shown for half the chain length (from i = 1 to 
32): the gap parameter is symmetric about the middle. 
In Fig. 1(a) we plot the result from the BdG equations, 
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while in Fig. 1(b) we show the corresponding results from 
the Anderson prescription. The first thing to note is that 
in either case the behaviour near the surface is markedly 
different as a function of electron density. At half-filling 
the gap parameter actually peaks at the surface, with 
several 'Friedel-like' oscillations ensuing towards the cen- 
ter of the sample, while at low fillings the gap parameter 
is much smoother by comparison. A comparison of the 
two figures shows quantitative differences, but overall, 
qualitatively they are very similar. It is evident that the 
Anderson prescription captures the essence of the BdG 
results remarkably well. 
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FIG. 2. Cross sections of Ai shown in Fig. 1 at density (a) 
n = 0.2 and (b) n = 1.0. The BdG and Anderson results 
are now plotted together, in solid and dashed lines, respec- 
tively. It can be clearly seen that the Anderson prescription 
reproduces the BdG results very well. 

The accuracy of the Anderson results can be seen more 
closely in Fig. 2, where the cross sections of Aj versus i 
in Fig. 1 for (a) n = 0.2 and (b) n = 1.0 are shown. It 
is clear that the essential features of the BdG results are 
reproduced in the Anderson approach. 

To examine the effect of impurities, we show in Fig. 3 
Ai as a function of i, for N = 32 with an impurity at the 
central site with varying energy (both negative and pos- 
itive). We have used PBC and an intermediate coupling 
strength, \U\ = 2t, at electron density n — 0.9. We have 
intentionally stayed away from half-filling, at which the 
ground state is not a superconducting state, but a charge 
density wave. One may recall that with periodic bound- 
ary conditions and with no impurities, the ground state 
for the attractive Hubbard model at half-filling is dou- 
bly degenerate: both superconducting and charge den- 
sity wave solutions coexist at this point. However, the 
presence of an impurity tilts the balance in favour of the 



charge density wave, and the BdG equations converge 
to a solution in which the pair potential, A,, is identi- 
cally zero at all sites. The Hartree potential, V i} on the 
other hand, oscillates as a function of site position. The 
Anderson prescription is unable to reproduce this (cor- 
rect) feature at half filling, and gives a superconducting 
solution with nonzero gap parameters. Also if the self- 
consistency of the Hartree potential is neglected in the 
BdG equations, this physics is missed, and the ensuing 
BdG result is similar to the Anderson solution. 
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FIG. 3. Gap parameter Ai as a function of site number 
i, for a chain of N = 32 with PBC and an impurity at site 
16 with varying potential energy, from the (a) BdG and (b) 
Anderson equations. Note that in (b) the scale of the gap is 
twice the scale in (b) . 

Returning to Fig. 3 we have plotted the BdG results 
in Fig. 3(a), while in Fig. 3(b) we show the results from 
the Anderson equations. The gap is suppressed at the 
site with a positive impurity potential (as is the site den- 
sity rii) and exhibits 'Friedel-like' oscillations around it. 
The Anderson prescription captures this behaviour qual- 
itatively, while it tends to underestimate the amplitudes 
of the oscillations (compare the solid curves in Fig. 3(a) 
and (b) for eiq — 0.5 1, and note the magnified scale in 
the latter graph). The Anderson results become better 
for stronger impurity potentials and for electron density 
n further away from half filling. For a negative impurity 
potential, when the potential strength is very weak, Aj 
(and rii) has a peak at the impurity site, as can be seen 
in Fig. 3(a) for ei6 = — 0.1 1 (the dotted curve). Though 
smaller in scale, the Anderson result in Fig. 3(b) has 
similar behaviour. As the potential strength increases, 
however, an attractive impurity tends to break the pair- 
ing, and suppresses the gap not only at the impurity site 
but also at surrounding sites (see the dashed curve for 
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ei6 = — 0.5 i in Fig. 3(a)). In such cases, the Anderson 
method overestimates the gap parameter around the im- 
purity site. This can be seen in Fig. 3(b), where the gap 
has the correct oscillating pattern, but with much smaller 
amplitudes. 
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FIG. 4. The Aj as a function of i, for a chain of A 7 = 32 
with PBC andeie = - 0.5t, for (a) \U\ = It and (a) |J7| = 3t 
The BdG and Anderson results are plotted in solid and dashed 
curves, respectively. 



good qualitative agreement with the BdG results. 

To our knowledge, the spatial dependence of the order 
parameter and the electron density has not been previ- 
ously explored in detail within the Anderson prescrip- 
tion. We have found, somewhat to our surprise, good 
agreement with results from the BdG formalism, for sur- 
faces and single impurities (aside from the weak scatter- 
ing limit). Throughout this study, it is important to note 
that in the vicinity of surfaces or impurities, "charge or- 
dered" states and superconductivity in general coexist 
(at the mean field level). In future work we will exam- 
ine various correlation functions and the local density of 
states, the latter of which has been (?J and will continue 
to be measured using scanning tunneling microscopy. 
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We study the attractive impurity case further in Fig. 4, 
where we show Aj and n, versus i for N = 32 with PBC 
and eie = -0.5t, for (a) \U\ = It and (b) \U\ = 3t. 
When the coupling is weaker than or comparable to the 
impurity potential, the Anderson approach captures the 
main features of the gap parameter around the impurity. 
This is the case for \U\ = It in Fig. 4(a), and indeed the 
Anderson results show excellent agreement with the BdG 
results. As the coupling becomes stronger compared to 
the impurity potential, the impact of the impurity be- 
comes more drastic, even with relatively weak strength. 
The gap is more suppressed around the impurity, and 
the density distribution exhibits 'Friedel-like' oscillations 
more enhanced. This can be seen for \U\ = 3 tin Fig. 4(b) 
(solid curves). On the other hand, the Anderson method 
yields a gap parameter that is more uniform as a func- 
tion of site position, and does not reproduce the correct 
oscillations in the site densities. 

In conclusion, we have formulated the BdG equations 
for a tight-binding model with an on-site attractive inter- 
action. We have retained the self-consistent Hartree po- 
tential in the BdG equations, and found that a single im- 
purity breaks the superconducting/charge density wave 
degeneracy which would otherwise exist at half-filling in 
this model. We have also formulated the prescription set 
out by Anderson, without impurity-averaging, and found 
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